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We review and interpret a modern approach to laser theory, steady-state ab initio laser theory 
(SALT), which treats lasing and amplification in a unified manner as a non-unitary scattering 
problem described by a non-linear scattering matrix. Within the semiclassical version of the theory 
the laser line has zero width as the lasing mode corresponds to the existence of an eigenvector of 
the S-matrix with diverging eigenvalue due to the occurrence of a pole of the scattering matrix on 
the real axis. In this approach the system is infinite from the outset and no distinction is made 
between cavity modes and modes of the universe; lasing modes exist both in the cavity and in the 
external region as solutions satisfying Sommerfeld radiation boundary conditions. We discuss how 
such solutions can be obtained by a limiting procedure in a finite box with damping according to the 
limiting absorption principle. When the electromagnetic and matter fields are treated as operators, 
quantum fluctuations enter the relevant correlation functions and a finite linewidth is obtained, via 
a generalization of SALT to include noise (N-SALT). N-SALT leads to an analytic formula for the 
linewidth that is more general than all previous corrected versions of the Schawlow-Townes formula, 
and can be evaluated simply from knowledge of the semiclassical SALT modes. We derive a simpler 
version of this formula which emphasizes that the noise is dominated by the fluctuations in the 
polarization of the gain medium and is controlled by the rate of spontaneous emission. 


I. INTRODUCTION 

The laser has been a critical enabler of the modern 
discipline of quantum optics via its highly monochro¬ 
matic and intense emission properties. The physical un¬ 
derstanding of steady-state laser emission has presented 
challenges since the early days after its invention, as 
first semiclassical QL1 and then fully quantum laser 
theories were developed M- Although these theories 
were able to capture most of the relevant physics, they 
tended to incorporate some phenomenological or approx¬ 
imate steps due to the complexity of dealing with the full 
non-linear and space-dependent semiclassical or quantum 
laser equations. Over forty years ago, several pioneers 
of early laser theory, Lang, Scully and Lamb, addressed 
some of the weak points in the theory in a classic paper 
entitled, “Why is the Laser Line so Narrow? A Theory of 
Single Quasi-mode Laser Operation” J'] • This work was 
aimed at a fundamental challenge in laser theory, dealing 
with the openness of the laser system, which invalidated 
standard Hermitian modal decompositions appropriate 
for closed systems. Prior to this work essentially all the¬ 
ories were framed in terms of “quasi-modes” of the type 
studied by Fox and Li and were based on phenomenologi¬ 
cal modal equations with cavity damping inserted to rep¬ 
resent the outcoupling of the modes Q . The goal of Lang 
et al. was to go beyond these quasi-mode equations for 
the laser field amplitudes to a more fundamental descrip¬ 
tion which included both the cavity and exterior regions 
explicitly (i.e. a model of “the universe”), and then work 
back from the modes of the universe to the quasi-modes 
of the cavity region alone. Using this approach the au¬ 
thors (1) “hope[d] to understand the mechanism leading 
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to the extreme monochromaticity of the laser radiation”, 
and (2) “to investigate the sense in which the present 
calculation, which does not include the cavity dissipa¬ 
tion as a phenomological loss mechanism, still leads to 
and implies a fluctuation-dissipation theorem” Q. 

While this approach to understanding the laser equa¬ 
tions was groundbreaking at the time, there has been 
much progress in laser theory in the subsequent years and 
it seemed to us a good time to re-examine these questions 
from a modern perspective. One reason for this is that 
during the past decade a new theoretical framework has 
been developed, known as the steady-state ab initio laser 
theory (SALT), which treats the openness of an arbitrary 
cavity exactly within the semiclassical theory , and 

thus has provided new insights into lasing phenomena. 
The resulting non-linear SALT equations for steady-state 
lasing can be solved efficiently using a non-Hermitian ba¬ 
sis set method and have been shown to agree to very high 
accuracy with brute force simulations of the atomic las¬ 
ing equations pjltHL including even N-level models with 
multiple independent lasing transitions pbj . The theory 
can also be generalized to treat injected signals (I-SALT), 
as will be relevant below [3- Moreover, recently Pick et 
al. have added Langevin noise to the SALT equations (N- 
SALT) and, using the SALT steady-state lasing modes 
as a starting point, calculated the quantum-limited laser 
linewidth analytically jT?}- Their analytic results have 
very recently been confirmed, again by direct simulations 
of the lasing equations, but with appropriate noise terms 
included [3 ■ It should be said that many aspects of the 
results obtained through SALT and N-SALT were antic¬ 
ipated by prior theoretical approaches, a partial list of 
which includes the work by Spencer and Lamb [3 - Fu 
and Haken [jtil - lzil . and Henry [3; however SALT and 
N-SALT appear to be the most general, flexible and in¬ 
clusive formalism yet developed for steady-state lasing. 



It treats the fields over the entire “universe” from the 
outset, but in a different manner from Lang et al ., and 
does not arrive at modal equations with damping in the 
conventional sense. We will now examine the questions 
raised by Lang et al. from the perspective of SALT and 
N-SALT. 


II. LASERS AS SCATTERING SYSTEMS 
A. Motivation and definition of the problem 

Most early lasers emitted from lasing inodes which 
were concentrated along some optical axis, e.g. the axis 
joining the centers of two mirrors or the geometric or in¬ 
dex profile defining a waveguide or fiber, and so there 
was a tendency to develop laser theories focused on one¬ 
dimensional or axially symmetric equations, and this is 
overwhel ming ly the type of formulation one finds in text 
books [24 h 26| |. However in the past two decades a great 
variety of laser systems have been developed with com¬ 
plex two or three-dimensional cavities, such as microdisks 
and microspheres [27l- l2j§ . micro-toroids [30], deformed 
(chaotic) disk cavities Him, photonic crystal defect- 
based and surface emitting lasers [341436^ and random 
lasers [37l - l43| . To describe such devices it was necessary 
to formulate laser theory in a manner that could both 
describe, and be implemented computationally, for cav¬ 
ities with an arbitrary geometry, and to include the full 
Maxwell wave equation, not its axial approximation. 

Moreover, some of the new lasers (e.g. random lasers) 
are extremely open systems, not even cavities in the 
usual sense; nor do they always emit well-behaved nar¬ 
row beams of light but rather generate complex emis¬ 
sion patterns in two and three dimensions. Thus it was 
crucial to find a description which did not over-simplify 
the spatial properties of the cavities and emission pat¬ 
terns, and which hence treated the cavity boundary and 
outcoupling exactly. The solution to this challenge pro¬ 
vided by SALT is to describe the laser as a scattering 
system. Any laser consists of two essential elements: a 
cavity region, given as a linear dielectric function, e c (x), 
which describes the scattering and confinement of pho¬ 
tons, either injected externally or internally generated, 
and a gain region, containing the medium which will be 
pumped to provide optical gain (sometimes cavity and 
gain medium essentially coincide). To include all cases 
the “cavity” will refer simply to a singly connected region 
large enough to enclose all the regions in which e c (x) ^ 1; 
and for simplicity we will assume the gain medium is fully 
contained within this region as depicted in Fig. [1] 

The laser so defined can function both as an amplifier 
of input radiation (when the gain medium is inverted) 
and as a self-organized oscillator/source of radiation in 
the absence of input, above the first threshold for lasing. 
From this point of view the laser is a scattering system 
which both scatters and amplifies input radiation, and 
also, at discrete frequencies, scatters specific internally 



FIG. 1. Schematic of lasers as scattering systems: a singly 
connected region surrounded by a last scattering surface con¬ 
taining an irregularly shaped cavity defined by a linear di¬ 
electric function e c ^ 1 and a non-linear gain region. Any 
incoming signals (blue arrows) are amplified and scattered by 
the cavity to outgoing states (red arrows); lasing corresponds 
to outgoing solutions without any input, and is only possible 
at discrete frequencies for a given cavity 


generated modes of outgoing radiation when sufficiently 
pumped. However unlike linear scattering systems, the 
laser is subject to non-linear saturation effects which are 
fundamental to its stability and must be taken into ac¬ 
count. 

Other non-linear effects within lasers, such as four- 
wave mixing, can destabilize their simple harmonic time- 
response, but these effects are neglected in SALT, and 
also (typically) in the earlier laser theories which moti¬ 
vated the work by Lang et al. 0 ■ Leaving out such effects 
is equivalent to assuming a stationary inversion density, 
giving rise to a gain susceptibility which does not mix 
frequencies. This assumption is always valid for single¬ 
mode lasing and amplification, which will be the focus 
of this work; it is also valid under certain rather gen¬ 
eral conditions for multimode lasing in microlasers, and 
SALT has been shown to apply quantitatively to such 
cases fl3til6| . 

Before proceeding, it should be mentioned that the in¬ 
put/output theory developed by Collett and Gardiner 
(and others) also considers the cavity to be a scatter¬ 
ing element [44], EH and is able to treat the full quan¬ 
tum operator properties of the scattered fields. This ap¬ 
proach has certainly contributed to our understanding of 
open quantum systems and the role played by the reser¬ 
voirs and noise processes. However, while there has been 
substantial effort expended to describe lasing within this 
framework [46j - [49l ]. it remains an open problem how to 
include the strong non-linearities present in active lasing 
systems in this formalism; whereas at the semiclassical 
level this is done quite accurately by SALT. 
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B. Scattering in the SALT formulation 

Initially we focus on the single-mode lasing case, as 
was done by Lang et al., so that the positive frequency 
component of the electric field will be given by: 

E (+ >(x,f) = ^(x)e" iwt , (1) 

where the real number, w, is either the input radiation 
frequency (for the amplifier) or the unknown frequency 
of the laser oscillator, which must be self-consistently de¬ 
termined. M'(x) is the unknown spatial lasing (or am¬ 
plified) mode, not assumed to be equal to any passive 
cavity mode or quasi-mode and is also self-consistently 
determined; such lasing/amplifier modes describe the so¬ 
lution over all space, not just inside the cavity. 

Assuming stationary inversion and neglecting quantum 
and thermal fluctuations, we can describe \l/(x) for an 
arbitrary laser by a non-linear wave equation of the form: 

[V x Vx + (e c + 4 t rx 9 )fc 2 ]¥(x) = 0, (2) 

where k = uj/c, and \g is the complex, frequency- 
dependent saturable gain susceptibility, which can be cal¬ 
culated for an arbitrary N-level atomic gain medium in 
steady-state [3] • For the case of a two-level gain medium 
it takes the form 

. , , = __ 
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(3) 

where Dq is the pump parameter (proportional to the 
inversion), qj_ is the dephasing rate of the polarization, 
oj a is the atomic resonance frequency, and T(u; — w a ) is a 
Lorentzian centered at u> a - Do and "E are measured here 
in appropriate dimensionless units [3 , and this form as¬ 
sumes that there is only one mode, either lasing or am¬ 
plified, active in the medium. (This is the only case we 
will consider explicitly here, but we will review the more 
general multimode SALT formalism in Appendix!©. The 
non-linear gain susceptibility represents the polarization 
of the gain medium and is the source of the amplifica¬ 
tion and (above threshold) of laser emission. The lack 
of any Langevin force terms in Eq. © is a result of the 
semiclassical approximation, in which the operators for 
the electromagnetic and polarization fields are replaced 
by their average values, eliminating the effects of quan¬ 
tum and thermal fluctuations. Such terms are taken into 
account in N-SALT, which will be discussed briefly in 
Sec. IVII1 but were neglected (intentionally) in the work 
of Lang et al. 

Assuming the system is pumped below the first lasing 
threshold and the inputs are small enough to generate 
negligible saturation of the gain, this equation defines a 
linear non-unitary scattering (S-)matrix for input fields. 
Depending upon the geometry there are different natural 
basis states describing input and output modes, Hi ^(x), 
Hi + *(x), from and to infinity, i.e. sufficiently far away 


from the cavity region, 

=E^ H L- ) ( x ) + ^ H i +) ( x )] 

Ot 

= 5> a (HM(x) + S a , ? Hj + >(x))] (5) 
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If the pump is strong enough to invert the medium and 
create gain at the input frequency, the S-matrix will be 
amplifying, with eigenvalues greater than unity. When 
functioning as an amplifier the boundary conditions are 
familiar: an input wave from infinity interacts with the 
cavity plus gain medium, scattering according to the ap¬ 
propriate dielectric or metallic boundary conditions in 
different regions of the cavity and giving rise to an out¬ 
going wave which differs from passive scattering in that 
it is amplified and more photon flux emerges from the 
scattering region than enters it. The case of large inputs 
leading to a saturated gain and the resulting non-linear 
S-matrix can be handled by I-SALT and will be included 
in the discussion below. 

As noted, In this formulation there is no distinction 
between cavity modes and modes of the universe, every 
amplified mode exists inside and outside the cavity as a 
steady-state scattering solution. However we now con¬ 
sider to the more challenging and subtle case of laser 
oscillation in the scattering picture. 


III. QUASI-MODES, RESONANCES AND 
QUASI-BOUND STATES 

Any scattering region of finite extent with a real po¬ 
tential or dielectric function, and hence with a unitary 
S-matrix, will have associated with it resonant solutions; 
i.e. stationary solutions at discrete complex frequencies 
which are purely outgoing. It is easiest to think of these 
as evolving out of the bound states of the system at real 
frequencies, Cj n , which would occur if the region were sur¬ 
rounded by a perfect lossless reflecting barrier. As the 
reflectivity of the fictitious barrier is reduced from unity 
to zero these real frequency solutions, corresponding to 
standing waves, evolve into complex frequency solutions 
with a substantial outgoing wave component (and no in¬ 
coming wave from infinity). Due to the purely outgoing 
boundary conditions it can be shown that the imaginary 
part of the frequency must be negative, and the S-matrix, 
continued to complex frequencies, will become infinite 
with simple poles at u> n — i') n /2, (y n > 0). Hence these 
purely outgoing solutions are discrete and countably infi¬ 
nite, in one to one correspondence with the bound states 
of the closed cavity. 

In the asymptotic region the frequency and wavevec- 
tor of the solutions are linked by the dispersion relation 
and the sign of will cause the corresponding spatial 
solution to grow exponentially at infinity and carry a di¬ 
verging outgoing photon flux. Hence these solutions are 
not physically realizable, but have physical significance 




because the S-matrix will have a resonant response near 
the real values to n , with a bandwidth given by y n (as¬ 
suming well-separated high finesse resonances). The cor¬ 
responding spatial solutions are referred to by various 
terms: quasi-modes, quasi-bound states and resonances. 
Historically the use of multiple terms seems to have arisen 
from the transition from 3D confinement, as e.g. in a mi¬ 
crowave resonator, where the connection to bound states 
is more apparent, to open, axial mirror-based resonators, 
where the resonances perpendicular to the mirror axis are 
essentially absent, and the Fox-Li method was developed 
to determine the high Q resonances. However there is 
no fundamental difference between the various types of 
outgoing solutions: they all correspond to poles of the 
passive cavity S-matrix. 

Such states also are very well-known from quantum 
scattering theory (where they are usually referred to as 
quasi-bound states), and are often considered to be phys¬ 
ical states because if an initial condition is given cor¬ 
responding to the quasi-bound state truncated to the 
“cavity region”, it will leak out to infinity with a rate 
given by y n (this property was explicitly demonstrated 
for the one-dimensional cavity treated by Lang et al.). 
However in this context we would like to emphasize that 
the true quasi-bound states, including the exponentially 
growing component at infinity, are not physically realiz¬ 
able steady-state solutions and cannot correspond to the 
lasing mode. Also, since they occur at the poles of the 
S-matrix for the passive cavity, they can only describe 
the properties of the passive cavity, independent of the 
gain medium and spectrum. Thus using “quasi-mode” 
also as a term for the lasing mode in the presence of gain 
can, and in our view often has, led to confusion. 
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FIG. 2. Simulations for a D-shaped laser cavity with R = 
rg = .5 R, n = 3.5 as shown in the bottom schematic, 
with a uniform two-level atomic gain medium characterized 
by \ a = 1 fim and yx = lOnm. (Top) Plot of the Lorentzian 
gain profile for the atomic medium. (Bottom) Plot of the 
motion of the poles of the S-matrix as the D-shaped cavity is 
pumped to the first lasing threshold. The pole corresponding 
to the first lasing mode is seen to reach the real axis (marked 
with a red circle), and the spatial profile of the lasing mode 
at threshold is shown just above. 



IV. THRESHOLD LASING MODES AND 
CONSTANT FLUX STATES 

However, as described above, the properties of the 
laser/amplifier are not governed by a unitary S-matrix 
once the complex gain susceptibility is taken into ac¬ 
count. This complex susceptibility violates the Hermitic- 
ity of the wave equation even for the closed cavity and 
invalidates the standard modal decomposition. But it 
makes possible purely outgoing solutions at real frequen¬ 
cies, which hence are physical solutions over all of space. 
The non-equilibrium steady-state gain susceptibility de¬ 
pends on the pump level, with its imaginary part chang¬ 
ing sign from positive (absorbing) to negative (amplify¬ 
ing) at the transparency point, at which the upper and 
lower lasing levels are equally occupied and the suscepti¬ 
bility vanishes. At this pump value the poles of the cavity 
with gain coincide again with those of the passive cav¬ 
ity. As the pump is increased beyond the transparency 
point, all of the passive cavity poles continue to exist, 
but typically they move up in the complex plane towards 
the real axis [01 [H(| HU , at a rate that depends on their 
position in the gain spectrum. The S-matrix of the am¬ 


plifier can be measured for real frequencies as the pump 
increases, and the height of the resonances correspond¬ 
ing to each pole increases (the transmitted flux is greater 
than the incident flux due to amplification), while their 
width decreases (gain narrowing), because the effective 
loss rate is proportional to the imaginary part of the pole 
frequency, which is decreasing. If there is no input to the 
amplifier then there is still no emitted field and no indica¬ 
tion within semiclassical theory of the approach to lasing 
(although when quantum fluctuations are included there 
will be enhanced noise, amplified spontaneous emission 
(ASE), near the gain center). If lasing is achieved, one 
of the poles of the S-matrix reaches the real axis, gain 
completely balances loss, and we have a zero width res¬ 
onance on the real axis. This process is demonstrated in 
Fig. [2] for a non-trivial two-dimensional cavity, a chaotic 
“D-shaped” dielectric cavity [33], where each pole is seen 
to move towards the real axis an amount proportional its 
distance from the center of the gain curve and one of the 
poles experiences enough gain to reach threshold. 

It first may seem that we are back at the familiar situ¬ 
ation of true bound states, which correspond to poles on 
the real axis, but this is not so. In the former case there 
is no photon gain and one simply has a decaying standing 








wave pattern, trapped in some region of space. Here we 
have a different type of object, which we call a “threshold 
lasing mode” (TLM). This type of solution never occurs 
in Hermitian theories, a purely outgoing solution at real 
frequency; it reflects the fundamentally non-Hermitian 
nature of the electromagnetic wave equations with a com¬ 
plex susceptibility, which is an effective theory in which 
the matter degrees of freedom are hidden. Such a solution 
is characterized by a steady-state photon flux out to in¬ 
finity, but with no unphysical growth of intensity outside 
the cavity and an emission pattern corresponding to the 
eigenvector of the S-matrix with a diverging eigenvalue. 
It is an example of a constant flux (CF) state: photons 
are generated within the cavity but are conserved outside 
the cavity, hence corresponding to a true physical lasing 
solution 9,11,12], In Appendix[A]we will show that each 
TLM is just one member of a complete biorthogonal ba¬ 
sis set of CF states which forms a natural basis set for 
describing the non-linear lasing solution above threshold. 

However we still have a conceptual problem: the lin¬ 
ear S-matrix for this amplifier is now infinite at the lasing 
frequency. Within linear theory any input would be in¬ 
finitely amplified, violating conservation of energy. The 
resolution of this dilemma is that infinitesimally above 
threshold the saturating non-linearity kicks in and stabi¬ 
lizes the solution. In fact the steady-state lasing am¬ 
plitude rises smoothly from zero and has no jump at 
threshold, just a discontinuous derivative. Even though 
formally the S-matrix has a pole fixed on the real axis 
above threshold, any input wave at the lasing frequency 
simply increases the gain saturation and doesn’t lead to a 
runaway of the output. Strictly speaking, any such input 
actually forces the pole off the real axis, as we will show 
with an example in Sec. Hence the non-linearity is es¬ 
sential for stablizing the linear CF solution. However, the 
actual solution near threshold is accurately described by 
the single TLM calculated from the linear wave equation 
with gain, since the saturation term is not strong enough 
to change the spatial form of the solution significantly. 

Returning to the fundamental question of modes of the 
universe vs. quasi-modes of the cavity, we summarize the 
SALT perspective. 

• If one treats lasers as unbounded scattering sys¬ 
tems from the outset, all solutions are defined both 
inside and outside the cavity region and there is no 
distinction between cavity modes and modes of the 
universe. 

• If the laser is functioning as an amplifier, the so¬ 
lutions are simply amplified scattering states de¬ 
termined by the boundary conditions in the cavity 
region and the degree of gain at the input frequency. 

• If the laser is functioning as a self-organized oscil¬ 
lator the solutions are purely outgoing states corre¬ 
sponding to an eigenvector of the linear scattering 
matrix with diverging eigenvalue due to the exis¬ 
tence of a pole on the real axis. Such solutions 


require the presence of the saturating non-linearity 
to be physically stable. 

• For any given cavity and gain susceptibility such 
purely outgoing states can only occur at discrete 
frequencies because they arise from moving the dis¬ 
crete poles of the passive cavity onto the real axis. 
Within semiclassical theory (neglecting quantum 
and thermal fluctuations of the fields) these solu¬ 
tions have precisely zero linewidth. 

We reserve the term “lasing mode” for these special dis¬ 
crete solutions. 

In summary then, our response to the question “why is 
the laser linewidth so narrow?” is that indeed this is due 
to the extreme limit of gain narrowing: only at discrete 
frequencies can the linear gain diverge and generate self- 
organized oscillation. As we will show in the next section, 
this zero-width resonance cannot be observed with an 
input signal, i.e. in amplification mode. However these 
perfectly monochromatic purely outgoing solutions can 
be found by SALT and have been verified by brute force 
FDTD simulations of the laser equations for many ex¬ 
amples [3 HIS- Such infinitely sharp lines are not of 
course obervable in true laser emission due to quantum 
fluctuations as we will discuss in Sec. 1VIII However the 
smallness of these quantum effects lead to the extraordi¬ 
nary narrowness of the resulting lines, often many orders 
of magnitude narrower than the passive cavity resonance 
linewidths. 

The effects of quantum fluctuations were intentionally 
not included in the work of Lang et al., who thus also 
found a delta function laser line within semiclassical the¬ 
ory. In their approach they find that this extreme nar¬ 
rowing of the passive cavity resonance can be regarded as 
a “mode-locking effect” in which all of the modes of the 
universe within the passive cavity resonance linewidth 
contribute to one lasing quasi-mode with a single fre¬ 
quency. This point of view emerges from a limiting pro¬ 
cedure in which a finite system with quantized electro¬ 
magnetic modes is studied and at the end the system 
size is taken to infinity in the presence of small fictitious 
absorption to damp out reflections from the boundaries 
as they tend to infinity. The conditions for this pro¬ 
cedure to reproduce the correct infinite system behav¬ 
ior are discussed in Sec. eh below. We note that Lang 
et al. expressed dissatisfaction with invoking gain nar¬ 
rowing to explain the laser linewidth, because in their 
view such arguments (at the time) did not take into ac¬ 
count the non-linearity of the problem. However the con¬ 
ceptual framework of SALT/I-SALT, outlined here, does 
take non-linearity into account, essentially exactly, for 
single-mode lasing. We will now illustrate the previous 
concepts with a concrete one-dimensional example, very 
similar to the one studied by Lang et al., in the process 
making all of our statements above mathematically pre¬ 
cise. 
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FIG. 3. Schematic showing the system studied by Lang et 
al. 31 which contains two regions, the cavity and the much 
larger “universe” (which also includes the cavity). Damping 
is added to the external region to prevent reflections from the 
boundary at x = —L, and as discussed below, this must be 
done by a particular limiting procedure in order to obtain the 
correct result for the laser. 


V. LASING MODES AS ZERO-WIDTH 
RESONANCES: A SIMPLE EXAMPLE 

The example studied by Lang et al. was a one¬ 
dimensional cavity of length L + l, with perfectly re¬ 
flecting mirrors at x = l,—L and a partially reflecting 
delta function mirror at x = 0, as shown in Fig. [3] The 
“cavity” is the region [0, l} and the “universe” is the re¬ 
gion, [— L, 0], where in treating lasing a small uniform 
absorption is added over all space, the limit L —> oo 
is taken, and lasing emission is in a single direction to¬ 
wards x = —oo. For such a system the passive cavity 
will have a reflection coefficient of modulus unity at all 
frequencies, and resonant reflection will be distinguished 
only by a long delay time in reflection and an associ¬ 
ated build-up of field intensity in the cavity. We prefer 
to make a small modification of the Lang et al. example 
and study a symmetric two sided partially transmitting 
(Fabry-Perot) cavity with laser emission to x = ±oo, 
and with passive cavity resonances distinguished by en¬ 
hanced transmission, reaching unity at the center, and 
with width, 7 C . As with the Lang et al. example, the 
electromagnetic part of the problem will be described by 
the scalar Helmoltz equation with gain. However, unlike 
Lang et al., no limiting procedure is applied; our system 
is infinite from the start and no fictitious absorption is 
added to avoid reflections from perfect mirrors far away 
from the cavity. Hence, consistent with our scattering ap¬ 
proach, we start from continuous solutions for the passive 
cavity + “universe”, and will impose outgoing boundary 
conditions at infinity on the lasing solutions. 

We begin by reviewing the properties of such a one- 
dimensional optical cavity. In Fig. [I] we plot the trans¬ 
mitted intensity of an incident signal through the passive 
cavity at frequency Wi„. As expected we observe a series 
of resonances of the cavity, spaced by the free spectral 
range, A = i t/cL, with a width given by the passive cav¬ 
ity decay rate y c = — c/(2 L) In (R 2 ), where L is the length 
of the cavity. The incident and reflected waves on the left 
side of the cavity satisfy a continuity boundary condition, 
with both incoming and outgoing components, whereas 
the transmitted wave leaving the right side of the cavity 
also is continuous but in addition is purely outgoing. 

Any purely outgoing monochromatic wave is described 
mathematically by the Sommerfeld radiation condition 
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FIG. 4. (Top) Plot of the transmitted intensity of an injected 
signal with frequency Win, into a passive, two-sided cavity with 
partially reflecting mirrors (R = 0.8) on both ends, as indi¬ 
cated in the schematic. (Bottom) Plot of the poles, w p of the 
scattering matrix of the same cavity, corresponding to purely 
outgoing solutions at complex frequency. The real freuency 
resonances of the cavity, where complete transmission is seen, 
are related to the poles of the cavity. 


lim \x\ d ^ (d\ x \-ik)ij)(x) = Q, (6) 
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in which x = x/\x\, w = ck, and d is the dimensionality 
of the system. In one dimension the external waves are 
simple plane waves (and not e.g. cylindrical or spherical 
waves) and the Sommerfeld condition can be imposed 
anywhere outside the boundary; thus at the boundary on 
the right side of the cavity the field must be an outgoing 
plane wave, satisfying, 

[{d x -ik)il>{x)\ x=L = 0. (7) 

As noted, lasing would correspond to a solution satisfy¬ 
ing a purely outgoing boundary condition on both sides of 
the cavity (or in all directions in higher dimension). If we 
attempt to impose this boundary condition on the pas¬ 
sive cavity we find that no solutions are possible at real 
frequencies, but there are discrete periodically spaced so¬ 
lutions at 

7r 

X>m = ttt + m A - *7 c/2 (8) 

ZcL 

as shown in Fig. 01 The discreteness of the frequencies is 
a direct consequence of the Sommerfeld condition which 




















can only be satisfied when there is constructive interfer¬ 
ence in the cavity. The fact that all the resonance fre¬ 
quencies have the same imaginary part is not generic, but 
is particular to one-dimensional cavities of this type. As 
noted above, the presence of these purely outgoing (pole) 
solutions leads to resonant transmission on the real axis 
centered at Re[u; m ] with width y c . Applying the Som- 
merfeld condition to the outgoing solutions e ±lkx , and 
the relation ui = ck implies that |\l/(a;)| 2 ~ e +7c ^ at 
±oo, the unphysical divergence we alluded to above. 

In order to find normalizable physical solutions which 
are purely outgoing (i.e. lasing solutions) we need to add 
the effect of the gain medium to the cavity dielectric func¬ 
tion in the form of \g as defined in Eq. ([3]) above. For 
this simple one dimensional scalar model (similar to Lang 
et al .), below threshold, the electric field satisfies 

[d 2 + {nk) 2 ]^(x) = 0, (9) 

with the Sommerfeld boundary condition, where now the 
index of refraction n = ^e c + An ry g = m — in 2 , and 
the gain non-linearity is neglected, because we are below 
threshold and these are not yet physically allowed solu¬ 
tions (they are simply modified active cavity resonances). 
Note that the solutions inside the cavity here take the 
form Acos(n:r) + Bsin(nx), i.e. they are trigonometric 
functions of a complex argument. This will remain true 
at threshold and is a key point which holds more gen¬ 
erally: the electric fields for the physical system with 
gain cannot be purely real; if they were real they would 
conserve photon flux and describe true standing waves, 
whereas lasing or amplified modes must have running 
wave components. Any model which finds real spatial 
solutions for lasing or amplified modes has made an ap¬ 
proximation. Such approximations have often been made 
in laser theory, including by Lang et al. @,M (see 
further discussion in Sec. ED. 

When we add the gain medium to the cavity and pump 
above the transparency point (Do > 0) the imaginary 
part of Xg corresponds to amplification. This has an im¬ 
mediate effect, causing the poles to move upwards, as 
discussed earlier and as shown for this specific cavity in 
Fig. [5] While the pole solutions are still not physically 
realizable, we can observe a change in the cavity prop¬ 
erties via the behavior of the transmission resonances at 
real frequencies. To calculate this effect accurately with a 
finite input wave we must now use the full saturable gain 
susceptibility of Eq. ©. The simulations shown here are 
performed using I-SALT which provides a nearly exact 
description of the steady-state single-frequency response 
of an active optical cavity to injected signals and any self¬ 
generated lasing signals [lg , including the spatial degrees 
of freedom. SALT and I-SALT are reviewed in more de¬ 
tail in Appendix 0 

As the pole rises towards the real axis the amplified 
resonance narrows. Due to the frequency-dependent sat¬ 
uration of the gain medium, as already noted, injected 
signals with different frequencies will induce somewhat 
different trajectories for the pole in order to maximize 



Frequency, w in (Top), Re[w p ] (Bottom) 


FIG. 5. (Top) Plot of the transmitted intensity of an injected 
signal as a function of frequency, Wm, illustrating line narrow¬ 
ing with increased gain. li n = Io = 0.01, and the two-sided 
cavity has partially reflecting mirrors (R = 0.8) at both ends, 
(see schematic). The pump on the cavity is increased from 
the transparency point, Do = 0 (red), to just above the first 
lasing threshold, Do = 0.0040 (blue), calculated using I-SALT 
fl6| . As can be seen, the width of the cavity resonance nar¬ 
rows during this process. The slight shift of the resonance is 
due to line pulling from central atomic transition frequency, 
cj a = 60, with yx = 2. At the transparency point this sys¬ 
tem is identical to that shown in Fig. [4] and the outermost red 
curve seen here is identical to that from Fig. [I] (Bottom) Plot 
of the motion of the pole of the associated cavity resonance, 
cj p , through the complex plane as the pump is increased. This 
is also calculated using I-SALT, and thus includes the satu¬ 
ration of the gain medium due to the injected signal, which 
prevents the pole from reaching the real axis even though 
the final pump value is above the first non-interacting modal 
threshold. For the specific realization of the motion of the 
pole shown, the injected field is taken at CJi n = 60.25, which 
corresponds to the peak in the blue curve. The pole is also 
seen to drift due to line pulling. Frequencies and rates are 
reported in units of c/L. 


the effective gain 0. These are essentially line-pulling 
or pushing effects of the gain medium, which for a sin¬ 
gle frequency are always towards the atomic transition 
frequency. In the active cavity with dispersion, the pre¬ 
cise width of the transmission resonance can no longer 
be defined as the distance of the pole from the real axis, 
as probing different frequencies will yield slightly differ¬ 
ent distances of the pole from the real axis. These line¬ 
pulling effects lead to the mild asymmetry in the series 
of curves shown in Fig. [5j 

The bluest (and narrowest) curve shown in Fig. [5] cor- 














Frequency, w in (Top), Re[w p ] (Bottom) 

FIG. 6. (Top) Plot of the transmitted intensity of an injected 
signal as a function of frequency Win, into the active, two- 
sided cavity with partially reflecting mirrors (R = 0.8) on 
both ends, as shown in Fig. [5] The pump on the cavity is 
held fixed just above the first lasing threshold, Do = 0.0040, 
while the incident intensity is decreased from /i n = Jo = 0.01 
(blue), to 7i n = 0 (red), where the cavity is operating as a 
free-running laser, calculated using I-SALT [jj3]. As can be 
seen, the width of the cavity resonance continues to narrow 
during this process, finally resulting in the lasing signal with 
infinitesimal width. Note, the outermost blue curve shown 
here is the same as that shown in Fig. [5] (Bottom) Plot of 
the motion of the pole of the associated cavity resonance, w p , 
through the complex plane as the injected signal is turned 
off for the specific realization of Win = 60.25. This is also 
calculated using I-SALT, and thus includes the saturation of 
the gain medium due to the injected signal, which prevents 
the pole from reaching the real axis until the injected signal 
is turned off. For different values of Wi n (not shown) the pole 
is able to reach the real axis before the injected signal is fully 
turned off. 


responds to a pump value just above the non-interacting 
first lasing threshold. This may seem initially puzzling; 
wasn’t the lasing threshold the point at which the pole 
reaches the real axis? That definition only applies when 
there is no input signal. With an input signal, even at 
the lasing frequency, the saturation of the gain medium 
prevents the pole from reaching the real axis, keeping 
the response finite, and the amplifed resonance retains a 
non-zero width. 

To see the true zero-width semiclassical lasing solu¬ 
tion we slowly turn off the incident signal with the pump 
held above the first lasing threshold. With this protocol 
we see that the transmission resonance width continues 
to narrow, as shown in Fig. [Gl while the pole continues 
to move towards the real axis. As the injected signal is 


turned off, the pole reaches the real axis and the cav¬ 
ity begins to undergo self-organized lasing action, shown 
as the reddest curve in the top panel of Fig. [G] This 
curve has precisely zero linewidth; no radiation comes 
out of the cavity except at this single frequency with¬ 
out input waves (assuming effectively zero temperature 
in the cavity and neglecting quantum fluctuations until 
Sec. EH. This effect can never be fully achieved through 
gain-narrowing of an amplified input signal, for the rea¬ 
sons we described above: gain saturation forces all poles 
off the real axis. Hence this effect does require correct 
treatment of the non-linearity as stated by Lang et al. 
However we have shown here that it is possible to treat 
the non-linearity exactly at each stage and explain the 
zero width of the laser line quantitatively. 

The model we have treated here is essentially the same 
as that treated by Lang et al. except for the pedagogi¬ 
cal difference of choosing a two-sided cavity in which one 
can study transmission resonances; in the next section 
we will study exactly the model of Lang et al. by a dif¬ 
ferent but equivalent method. As noted, in the standard 
SALT approach the system is infinite from the outset and 
there is no dissipation outside the cavity, which is in our 
view the most natural model for a real laser system. In 
contrast, Lang et al. place their cavity in a much larger 
box which eventually will be taken to infinite size, and 
add dissipation external to the cavity, intended to damp 
unphysical reflections from the boundaries of the larger 
fictional box. We will discuss the effects of this difference 
in the next section. 


VI. LASING MODES AND THE LIMITING 
ABSORPTION PRINCIPLE. 


While it is clear that lasing modes should satisfy the 
purely outgoing (Sommerfeld) boundary condition, this 
leads to a rather unfamiliar description of the spatial be¬ 
havior of the lasing modes in terms of functions that are 
not power orthogonal, but instead satisfy a biorthogo¬ 


nality condition 25|j_ 56|, |57(. Basis functions of this type, 
the CF states Sua , are central to the standard com¬ 
putational algorithm of SALT. The approach taken by 
Lang et al. of beginning with a finite system and adding 
damping to remove the effect of the boundaries of the 
system is an alternative which has evolved into the well- 
developed computational technique of surrounding a cav¬ 
ity or region of emitters with a perfectly-matched (damp¬ 
ing) layer (PML). When done correctly this can repro¬ 
duce the imposition of Sommerfeld conditions to arbi¬ 
trary accuracy. In fact, newer computational methods 
for solving the SALT equations [5f| make use of PMLs 
and do not impose the Sommerfeld boundary conditions 
in the straightforward manner of matching to outgoing 
solutions at a last scattering surface as did the earlier 
methods of satisfying the SALT equations The 

PML approach gives essentially perfect agreement with 
the Sommerfeld approach to SALT. 














However, as noted above, the correct lasing solutions 
must be complex functions in space, so as to represent 
the non-zero photon flux, which grows inside the cavity 
and remains constant outside. For the simple ID cav¬ 
ities studied here and by Lang et al. , this means that 
the spatial solutions inside a constant index cavity must 
be trigonometric functions of complex argument, cor¬ 
responding to exponential growth towards the mirrors 
as well as oscillation. However in the work of Lang 
et al. they find a lasing mode of the form E(x,t) ss 
Coe ( '°‘~‘ lu ' ,t sin[fco(a; — l )] for 0 < x < l, where a is the dif¬ 
ference between gain and loss, ui is the lasing frequency, 
Co is a constant, and ko = Cl/c is the real wavevector cor¬ 
responding to the passive cavity resonance frequency, Cl. 
Since the spatial part of this function is real, it does not 
capture the photon flux through the partially transmit¬ 
ting mirror and is only an approximation to the complex 
sine wave solution which is correct for this one-sided cav¬ 
ity. 

In our view the origin of this problem is the neglect of 
the Limiting Absorption Principle [59l - [6l| , which gives 
the correct mathematical procedure for recovering the 
Sommerfeld radiation condition from a finite system. 
Lang et al. solved the spatial part of the problem for finite 
L. finding the quantized Hermitian basis states for the 
whole passive system (“modes of the universe”) and then 
addressing the lasing problem within this Hermitian basis 
set. Effectively the spatial part of the problem is fixed 
at finite L , and a small absorption, 7 , is subsequently 
added everywhere in the universe so as to eliminate the 
effect of reflection from the fictitious perfect mirror at 
x = — L. For this procedure to work one needs to solve 
the spatial problem with the finite damping, and carefully 
take the limit L —> 00 prior to 7 —> 0; this requirement 
is known as the Limiting Absorption Principle (LAP). 
From a practical computational point of view the LAP 
means that one needs to solve the non-Hermitian spa¬ 
tial problem with Diric.hlet (or other linear homogeous) 
boundary conditions and finite damping, while maintain¬ 
ing the condition 7 kL 1. The solutions so obtained 
have only exponentially small sensitivity to the presence 
of the boundary and reproduce to high accuracy those 
obtained by actually imposing directly the Sommerfeld 
conditions. Hence they are significantly different from 
the Hermitian solutions calculated before the damping is 
added (or when the limits are taken in the opposite or¬ 
der), which can carry no net flux as they are bound states. 
The modern PML techniques are a means for introduc¬ 
ing non-uniform absorption which is impedance matched 
so that reflections are negligible even when L is not very 
large and one reproduces the Sommerfeld solutions with 
a much smaller and manageable computational cell. But 
the principle is basically the same: in order to get the cor¬ 
rect spatial solutions one must treat the non-Hermitian 
problem with damping from the outset. 

To illustrate this explicitly we briefly analyze the Lang 
et al. system from this point of view. The system is the 
one-sided cavity between [ 0 , l] attached to the “universe” 


between [— L, 0], as shown in Fig. [3j with perfect mir¬ 
rors (Dirichlet boundary conditions), E(l) = E(—L ) = 0, 
where E(x,t) is the electric field of the system. Before 
adding the gain and considering lasing, we explicitly in¬ 
clude the weakly absorbing material uniformly in the uni¬ 
verse, and focus on the region external to the cavity. The 
wave equation in this region is written as, 


a? - i± ±L d ; 


E(x,t) = 0, — L < x < 0. (10) 


We expand the field inside of the system in terms 
of a set of basis functions, Uk(x), so that E(x,t) = 
U k {x)e ~ iut , and solve for the form of these basis func¬ 
tions in the universe region as, 


U k (x) = c hk e ikx ~l kx + c 2ik e~ ikx+ * kx , 

L<x< 0, (11) 

where we have expanded the exponent for 7 C 1. The 
coefficients c\ yk and c 2tk and the quantized values of k 
are set by the boundary conditions U k {—L ), U k {l) = 0, 
and the matching conditions at x = 0 , and will vary with 
L 1 l and the transparency of the mirror at the origin. 

We first apply the LAP to obtain the Sommerfeld so¬ 
lution for this problem, taking the limit L —> 00 prior to 
7 —> 0. Performing the limits in this order on the solution 
of U k (x), Eq. (fill) , one finds that to satisfy the Dirich¬ 
let boundary condition, U k (—L) = 0, we must choose 
ci,fc = 0 , as the term multiplying this coefficient is di¬ 
verging as — L —> —00 for finite 7 . Then, with the spatial 
limit taken, we can safely reduce the absorption in the 
problem to zero, recovering 

U k (x) = c 2 , k e~ lkx , -00 < x < 0, (12) 

finding that all of the solutions of our problem are purely 
outgoing from the cavity, consistent with the Sommerfeld 
radiation condition, Eq. ©. However, if these two limits 
are taken in the opposite order, setting the absorption to 
zero before letting the length increase to infinity, we find 
that 


U k (x) = b k sm(k(x + L)), —L<x<0, (13) 

which is exactly the result of Lang et al. (outside the cav¬ 
ity region). With this result the matching to the cavity 
allows for real sine wave solutions, whereas with the out¬ 
going solutions only complex sine waves are allowed in 
the cavity. The complex sines provide the outgoing flux 
at the partially transmitting mirror which is needed for 
continuity with the field outside the cavity. 

These conclusions can be further reinforced by directly 
solving for the eigenfunctions of the finite system stud¬ 
ied by Lang et al., shown in the three panels of Fig. [7] 
As L is increased, the modes of this finite system sepa¬ 
rate into two sets. One set consists of discrete isolated 
solutions with smaller imaginary frequency (damping); 
the number of such solutions does not change with in¬ 
creasing L. The second set has larger damping, and gets 





denser and denser, forming a psuedo-continuum as L is 
increased. The first set corresponds precisely to the com¬ 
plex resonance frequencies of the passive cavity, as can be 
confirmed by plotting the known resonances of the cavity 
given by Eq. ([ 8 ]), shown in Fig. [Tfc). The spatial func¬ 
tions corresponding to the passive resonances are identi¬ 
cal to the resonant solutions of the infinite system with¬ 
out damping inside the cavity , as there is no loss there, 
and the boundary condition at the mirror is the same. 
Outside the mirror these solutions experience a damp¬ 
ing not present in the infinite system and hence decay 
slowly as they propagate outwards. On the scale shown 
in the top panel of Fig. [ 8 ] one sees only that the out¬ 
going boundary condition is satisfield and not the weak 
decay. In contrast, the spatial solutions corresponding 
to the dense set of damped modes of the universe can 
be seen in the bottom panel of Fig. [ 8 ] they form stand¬ 
ing waves outside the cavity, with small leakage into the 
cavity, and do not satisfy the purely outgoing conditions. 
Hence simply discarding the latter and keeping the first 
set is equivalent to finding the Sommerfeld solutions of 
the passive problem without damping, in the vicinity of 
the cavity. 

Finally, we can make contact with the SALT approach 
by adding gain (taken to be uniform for convenience) in 
the cavity region while leaving the damping present out¬ 
side the cavity. The added gain has a much bigger effect 
on the cavity resonance solutions which have high ampli¬ 
tude there, than on the external standing wave solutions. 
For large enough gain one of the cavity resonances can 
reach the real axis and will then correspond to a threshold 
lasing mode for uniform gain, as shown in Fig. [TJ; (green 
crosses). This is the essence of the method adopted in 
Ref. [58j except in that work a PML damping layer is 
used for higher computational efficiency and the saturat¬ 
ing non-linearity is taken into account. Thus the Lang 
et al. method agrees with the SALT approach if care is 
taken to apply the limits in the right order. 


VII. QUANTUM LIMITED LASER 
LINEWIDTH: AN ALTERNATIVE FORM 

Since this paper is devoted to the “SALT perspective” 
on the laser linewidth, we include here a section on the 
finite laser linewidth due to quantum fluctuations, which 
has recently been addressed within two different gener¬ 
alizations of SALT, one based on input-output theory 
nam and another based on non-linear perturbation 
theory around the SALT solution E3- The fundamental 
origin of a finite quantum linewidth for the laser is the 
presence of spontaneous emission, which leads to ran¬ 
dom jumps in the time-dependent phase of the output 
field. The resulting phase diffusion is linear in time and 
the phase diffusion coefficient determines the intrinsic 
linewidth. The basic physical understanding of this goes 
back to the seminal work of Schawlow-Townes f(32| , which 
was subsequently substantially extended by Lax, Henry, 
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FIG. 7. Plot of the resonances of the finite system with exter¬ 
nal damping depicted in Fig. [3] for three different values of L, 
L = 25 1, L = 100/, and L = 200/, from top to bottom, except 
that the optical confinement is instead provided by a passive 
cavity dielectric of n = 6 in the cavity, rather than a par¬ 
tially reflecting mirror. In the external region an absorption 
of 7 = 0.003 is present. In (c) the isolated solutions of the 
problem with damping are compared to the resonances of the 
same cavity subject to the Sommerfeld boundary condition 
on an infinite system without damping, given by Eq. (JH|) (red 
circles). Perfect agreement is found as discussed in the text. 
(Note that due to the presence of the passive cavity dielectric 
constant the right side of Eq. (J 8 |) must be multiplied by 1/n in 
making the comparison). The green crosses shown in (c) show 
the locations of the resonances when uniform gain is added to 
the cavity to find the threshold lasing mode, n = 6 — .005i. 
As can be seen, the poles corresponding to resonances of the 
cavity have moved towards the real axis,with the rightmost 
one reaching the lasing threshold, while those corresponding 
to the resonances of the universe have stayed relatively fixed, 
as discussed in the text. This is thus an alternative method 
to implement the Sommerfeld conditions for lasing modes. 



■*— universe cavity —*■ 


FIG. 8 . Plot of the modal amplitudes for a small region near 
the cavity-universe boundary for two modes, one correspond¬ 
ing to a resonance of the cavity (top), and the other corre¬ 
sponding to a resonance of the external universe (bottom). As 
can be seen, the cavity mode is a nearly pure running wave 
after it exits the cavity. 













Petermann and many others. 

Traditionally, the Schawlow-Townes (ST) linewidth is 
written in terms of the cavity decay rate, j c , as 


Slost = 


Tm oil 
2 P Hsp ' 


(14) 


where coo is the lasing frequency and P is the total out¬ 
put power. Here we have chosen to include in the ST 
formula the “incomplete” inversion factor, n sp , which is 
sometimes treated as a correction to it (n sp > 1 mea¬ 
sures the degree of inversion of the laser frequency and 
tends to unity when the lasing transition is completely 
inverted, the limit assumed by Schawlow and Townes in 
their original work) [63j . 

As discussed extensively above, semiclassical laser the¬ 
ory neglects quantum effects such as spontaneous emis¬ 
sion and leads to a delta function linewidth. In order to 
go beyond this one needs to treat all three fields, the elec¬ 
tromagnetic field, the polarization, and the inversion as 
quantum operators and consider their higher moments. 
The laser linewidth will be obtained from the correlation 
function of the electric field near the lasing frequency. It 
is well known that the effects of these quantum fluctu¬ 
ations can be captured by adding appropriate Langevin 
(random) noise terms to the equations of motion for the 
three fields [§4j ■ These noise terms will be related to the 
dissipation and amplification in the system by the appro¬ 
priate fluctuation-dissipation relations. This was one of 
the questions raised in the original work by Lang et al. It 
is now understood that semiclassical theory violates these 
relations but they are restored by including the Langevin 
terms. The recent work by Pick et al. [17| provides the 
most complete treatment of these noise terms and the 
formalism they develop is referrred to as N-SALT (SALT 
+ noise). 

Pick et al. begin with the time-dependent wave equa¬ 
tion for the field and add the appropriate noise term 


[V x V x E(x, t) + e c E(x, t)] = 47 rP + F tot , (15) 


However, in this section we would like to make a differ¬ 
ent conceptual point, not discussed in the earlier papers, 
relating to the interpretation of the Schawlow-Townes 
formula. It seems initially puzzling that the noise due 
to spontaneous emission in the gain medium should lead 
to the standard ST formula, which is expressed only in 
terms of the passive cavity linewidth, 7 c , the total out¬ 
put power, and the lasing frequency. One would expect 
the spontaneous emission rate to appear as a fundamen¬ 
tal quantity, not the cavity decay rate; and here we will 
show that the ST linewidth can (and perhaps should) be 
written in an equivalent form where this is the case. 

The total Langevin force, F tot {t ), in units of 
th e qua ntized field amplitude operator, F tot = 
i\/2irtvjjF to t (t) Vi' 0 (x), has a direct contribution from the 
amplitude fluctuations of the field, F a (t), due to quan¬ 
tum fluctuations from outside the cavity, but has also 
has indirect contributions from Langevin forces driving 
the atomic polarization and inversion, F p and F r i respec¬ 
tively [64j , 


Ftotit ) 


F a {t) + 


1 


loq- u} a + ill. 


'y ] 9aF p (t) 


(16) 


where the summation runs over all of the atoms in the 
cavity at locations x a and 


g a = id- $o(x a )^/^, (17) 

is the dipole coupling element for the quantized elec¬ 
tric field. 6 is the dipole matrix element for the las¬ 
ing transition, and tFo^) is the normalized mode profile, 
mi 2 dx = 1, of the electric field at frequency loq. The 
correlation function for the total effective Langevin force 
can then be calculated Q as, 


(Fl t (t)F tot (t'))=S(t-t') 


[lenth 


in which F fot (x, t) is the total effective noise driving the 
field in the cavity in units of electric field. Since this ran¬ 
dom drive is broadband, it is no longer possible to assume 
a monochromatic steady-state exists as we do in semiclas¬ 
sical theory. However only the noise near the lasing fre¬ 
quency is strongly amplified and by expanding the field 
around the SALT steady-state solution one can derive 
non-linear coupled mode type equations which permit ab 
initio calculation of the relevant correlation functions and 
hence of the most general form of the laser linewidth yet 
derived B3- This form includes all of the important cor¬ 
rections to the ST formula: the correction due to relax¬ 
ation oscillations (a factor) (23 l 65L due to openness of 
the cavity (Petermann) factor]56, 57,®il , due to dis¬ 
persion (bad-cavity factor) [631 7 , and corrections 

not previously studied due to the spatial variation of n sp 
[63[ . We will not describe this formalism in any more 
detail here, but will briefly reproduce the essential result 
at the end of this section. 


E 


1 9o 


7i 


7±(1 + <da» + -^-(do ~ {da )) 


(18) 


in which nth is the number of thermal photons at the 
lasing frequency, 7 _l is the dephasing rate of the polar¬ 
ization, and d a is the inversion of the ath atom, which 
would relax to the steady-state value of d 0 in the absence 
of an electric field at the rate 7 ||. We have also simplified 
by assuming the lasing signal to be on resonance with the 
atomic transition, loq = co a , and that the dipole matrix 
element and field are co-linear, 9 ■ \Fo = fl’I'o- 

The first term in Eq. (TTSl) corresponds to the contribu¬ 
tion from thermal and zero point fluctuations from out¬ 
side the cavity to the field inside the cavity; as such it 
is controlled by the decay rate of the cavity, i c . The 
second and third terms stem from spontaneous emission 
events leading to the dephasing of the atomic polariza¬ 
tion, and these terms depend on the parameters of the 
gain medium. For optical frequencies, the number of 






















thermal photons present is negligible, so the direct field 
fluctuations do not have any effect on resul ting linewidth. 
Moreover, for most lasers (class A and B) [73}, 7 y <C 7 ±, 
so that only the second term in Eq. ca remains, 


As noted above, the coefficient of the Langevin noise 
determines the laser linewidth by standard manipula¬ 
tions. Using the more correct form of this coefficient 
above leads to a quantum-limited linewidth, 


(FUt)F tot (t'))=5(t-t / )—J2\9 a \ 2 (l + ( d a)). (19) 


$ _ I lsp(%)D(x)n sp (x)dx 

2 rip 


(25) 


Assuming that the atoms are uniformly distributed over 
the cavity, the sum over the atoms in the cavity can be 
converted into an integral, 


N 2 (x) 

N 2 (x)-N 1 (x) 


f 

7-L 


D(x)\'i’o(x)\ 2 dx, (20) 


where N 2 (x) (Ni(x)) is the density of atoms in the upper 
(lower) lasing level, and g a = g , '^Q(x a ). Equation (12(1 
now expresses the Langevin force in terms of the incom¬ 
plete inversion factor 1631 ]. 


°sp 


(x) = 


N 2 {x) 


n 2 (x)-N! ( x y 


( 21 ) 


Traditionally at this step the atomic populations are 
approximated to be uniform in space which clamps the 
inversion to its value at the lasing threshold, D(x) ss 
Dq = DQ hr ' > j34|. In this approximation, the integral in 

Eq. (l20l) becomes just ^n sp , and can be related back 

to the cavity decay rate by 


D 


(thr) 
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Iff? 


7-L 


7c 
2 ’ 


( 22 ) 


which states that at the lasing threshold, gain must equal 
loss. Using this approximation contributes one factor of 
7 C in the standard ST linewidth formula, Eq. m • How¬ 
ever this is a rather crude approximation. Due to spa¬ 
tial hole-burning the atomic populations generally vary 
substantially in space, and as a consequence the overall 
inversion is not clamped at the threshold value. 

Using Eq. (El obscures the fact that the physical 
origin of the noise is spontaneous emission in the gain 
medium. It then is more natural to relate the linewidth 
to the rate of spontaneous emission , 7 sp , which can be 
calculated from the rate of stimulated emission, 7 se m. 
as 

/ \ 7 SE^x') 2 |g | . / \|2 

7 sp{ x ) — ^ — r^o?) 5 (23) 

2 n p 7± 


where n p is the number of photons in the steady-state las¬ 
ing field, and is related to the intensity by Iq = 2TihuJon p . 
Thus, we can rewrite Eq. El in terms of the rate of 
spontaneous emission, keeping the full space dependence 
of the problem, as 


{F} ot {t)F tot (t')) = 5(t-t') 


7 sp (x)D(x)n sp {x)dx . (24) 


Typically, the next step is to relate this to the total out¬ 
put power of the cavity via 

Pst = hlOJcUp, (26) 

which is also only an approximate relation, neglecting 
spatial variations of the intensity. Using both of the ap¬ 
proximations mentioned above and some other simplifi¬ 
cations valid for high Q cavities (and neglecting the am¬ 
plitude contribution contained in the a factor), one then 
arrives at the standard ST linewidth form, Eq. El. 

However, the actual power emitted can be calculated 
using the true spatially-varying fields using Poynting’s 
theorem [ 7 H, as 

-P= ^^/ Im [-e]|^°(a;)| 2 dx. (27) 

Using specific form of the imaginary part of the dielec¬ 
tric function for a two-level gain medium (and using real 
units, not SALT units), the linewidth can then be written 
as 


Suj = 


hjj 0 

2 P 


’y S p{x)D(x)n sp (x)dx / ■y sp (x)D(x)dx. 


( 28) 

Here we have derived a more physically intuitive and 
accurate form of the laser linewidth formula in which the 
critical role of spontaneous emission in the gain medium 
leaps out: the linewidth is determined by the spatially- 
varying, mode-dependent spontaneous emission rate and 
by the steady-state saturated inversion in the presence of 
spatial hole-burning. Only if one makes the approxima¬ 
tion of f "/spDdx « 7 C does this reduce to the standard 
Schawlow-Townes formula, and in making this substitu¬ 
tion one tends to hide the true origin of the linewidth, 
which is spontaneous emission. 

While the above argument captures the basic physics of 
the ST formula in an improved form, it is not a rigorous 
general result, as is the formula derived by Pick et al. 

17], which includes essentially all known corrections in 
a more accurate form. However the N-SALT linewidth 
does have the same general features as our simpler result, 
Eq. EH]), and for a single lasing mode can be written as 


r &*>() _ 2 

OWN-SALT = 77^7(1 + a ) 


w o / Im ? I^opdx f Im[e]^^|'4'o| 2 dx 


(29) 












where S&q is the normalized SALT lasing mode, and a 
is a generalization of the Lax-Henry a factor. Very re¬ 
cently this formula has been quantitatively confirmed us¬ 
ing first-principles numerical simulations [18[ for laser pa¬ 
rameters where it predicts significant departures from the 
standard (corrected) form of the Schawlow-Townes for¬ 
mula. Note that in Eq. (l29l) the cavity decay rate is com¬ 
pletely absent, while the imaginary part of the response 
of the gain medium plays a critical role. 

VIII. CONCLUSION 

In summary we have presented a scattering framework 
for semiclassical laser theory and showed how precisely it 
leads to zero laser linewidth, due to the discreteness of the 
poles of the scattering matrix. This approach treats the 
full system as infinite from the beginning and does not 
introduce a finite “universe” and impose a limiting pro¬ 
cedure and damping to recover the infinite system result, 
although such an approach is a valid alternative method 
if the limiting absorption principle is observed. We were 
inspired to analyze this question by the classic work of 
Lang et al. which introduced the notion of “modes of the 
universe”. Our scattering point of view underlies recently 
developed ab initio laser theory (SALT), which leads to 
efficient computational tools for calculating steady-state 
lasing properties. A generalization of the ab initio theory 
to treat noise allows accurate calculation of the intrinsic 
linewidth due to quantum fluctuations when one goes be¬ 
yond semiclassical laser theory. We emphasize that the 
quantum-limited linewidth depends directly on the prop¬ 
erties of the gain medium in the active lasing state, and it 
may be useful to write it in a form which does not involve 
the passive cavity decay rate, but rather emphasizes the 
spontaneous emission rate. 

ACKNOWLEDGEMENTS 

We acknowledge partial support by NSF grant DMR- 
1307632. ADS thanks Marian Scully for drawing our 
attention to the work of Lang et al. We thank Steven 
Johnson and Adi Pick for very useful conversations about 
the laser linewidth, noise effects, Perfectly Match Layers 
and the limiting absorption principle. We would also like 
to thank Li Ge for stimulating discussions on methods 
for finding the poles of the S-matrix for two-dimensional 
cavities. 


Appendix A: Review of SALT 

In this Appendix we will provide a brief review of SALT 
to assist in understanding the discussion in Sec. hyi and 
the results presented in Sec. However, the interested 
reader is referred to works dedicated to the SALT and 
the CF basis nnnm , and its elaborations to treat 


amplified modes [l(|, and complex gain media mm. 
The SALT equations find the steady-state solutions to 
the semiclassical Maxwell-Bloch equations, for which the 
wave equation was given above in Eq. ©. For this brief 
review we will only consider the simplest case of two-level 
atomic gain media, which leads to the Bloch equations 
for the evolution of the atomic polarization: inversion are 

id 

5 f P+(x, t) = - ( 7± + iu a ) P+ - - (0 • E) 6* , (Al) 

Wx,t) = - 7|| (D - D 0 (x)) - | (P" - P+) • E. 

(A2) 

For an electric field with a single frequency component, 
or operating in the stable multi-mode regime, we can 
expand the electric field over its N constituent frequency 
components as 

N 

E+(x,f)=^E CT (x)e- i ^ t , (A3) 

a 

where E + is the positive frequency component of the elec¬ 
tric field, E = 2Re[E + ]. In this expansion we make no 
distinction between lasing inodes, which are associated 
with the poles of the scattering matrix, and amplified 
modes due to injected signals, which are externally gen¬ 
erated. 

The critical assumption in SALT is the stationary in¬ 
version approximation (SIA), which states that dtD = 0, 
even when multiple frequencies in the electric field are 
present; i.e. the populations do not respond to beat¬ 
ing terms between the modes, that would lead to four 
wave mixing and frequency comb generation. This as¬ 
sumption requires that the non-radiative relaxation time 
of the gain medium is much slower than these beating 
terms, 7 || <C A, where A is the minimum frequency 
difference A = present in the set of frequen¬ 

cies comprising the electric field, {w CT }, and is satisfied 
by nearly all microcavity lasers. In the parlance devel¬ 
oped for the study of temporal chaotic laser behavior, 
the_SIA is valid for class A and B lasers, but not class 
c [zl|7l. For class B lasers, which include most semi¬ 
conductor based devices, the mode spacing must also be 
well separated from the relaxation oscillation frequency, 
so as to avoid resonantly driving relaxation fluctuations 
which could destabilize the multimode solution. The re¬ 
laxation oscillation frequency is wro ~ ^/k7_l [z 3 , where 
k = 7 c/2 is the field decay rate of the cavity. For mi¬ 
crocavities without an injected signal, k < Afsr, the 
mode spacing from the free spectral range of the cavity, 
so wr,o < y / Afsr7|| < Afsr, and as such, the SIA is 
still valid when 7 || < Afsr jlSj- Furthermore, SALT 
has been rigorously tested using direct finite-difference 
time-domain (FDTD) simulations of the Maxwell-Bloch 
equations, with excellent quantitative agreement found 

EH®. 

As the SIA removes any time dependence from the 
inversion, the polarization must now contain exactly the 




same frequency components as the electric field, 

N 

P+(x, i ) = ^ P(T (x)e-^. (A4) 

G 

Using the multimode expansion for the electric field, 
Eq. (IA3I) . the polarization equation, Eq. <EH>. can be 
solved independently for each frequency component, pro¬ 
viding an expression for in terms of D and and 
inserted into the wave equation as 


V x V x - e c (x) + 


47 r| 0 | 2 Z)(x) 
h(cOcr - UJ a + *7_1_) 


K E ct (x) 
= 0, (A5) 


which explicitly identifies the gain due to the inversion of 
the atomic medium, and in which 6 has been taken to be 
collinear to the electric field. Likewise, using Eq. (1A2[) , 
we can solve directly for the inversion as 


D(x) = 


-Do(x) 


1 + 


4|9| 2 

fi 2 7_L7|l 


Ef r » 


|E„ 


(A 6 ) 


in which the denominator contains the effects of gain- 
saturation due to the electric field within the cavity, 
and r fj = 7 — ui a ) 2 + 'yjj is the Lorentzian gain 
curve. Together, Eqs. (IA5I) and (IA 6 I) form the fundamen¬ 
tal SALT equations (for the two level gain case): N self- 
consistent, differential equations for the unknown spatial 
profiles of the N modes in the electric field at N unknown 
frequencies, coupled through the non-linear spatial hole¬ 
burning in the inversion. These are to be solved with out¬ 
going Sommerfeld boundary conditions at infinity, which 
can only be satisfied at discrete frequencies. The SALT 
equations also reveal the natural units for the electric 
field and atomic inversion, -Esalt = 2|0|/(/t v /7_i_7||)-E’, 
and Dsalt = (4n\6\ 2 /hy±)D. 

The numerical solution of the SALT equations is per¬ 
formed by expanding the mode profile functions over a 
basis, 


EUM = X] a n <7)f "( X;a;<7 )> ( A7 ) 

n 

and then using a non-linear solver to find the expansion 
coefficients a! r f ' 1 . There are two different basis sets that 
have been used to solve the SALT equations. The first, 
and more developed, is based on a set of constant flux 
(CF) states @, [HI, which are defined as 

[-V x V x + (e c (x) + ??„F(x)) k 2 ] u n (x; uj) = 0, (A 8 ) 

and satisfy the outgoing boundary condition at frequency 
w at the surface of last scattering. Here, rj n (u}) and 
u ra (x;w) are the eigenvalue and eigenstates of the CF 
equation, and D 0 (x) = D 0 F(x) is the profile of the 
pump. A key aspect of the SALT approach is that for 
any given frequency (wavevector) the CF states form a 


complete biorthogonal basis for an arbitrary purely out¬ 
going solution; hence they can expand any solution to 
the SALT equations. 

Moreover, by observing the similarities between 
Eq. (IA 8 I) and Eq. (1A5I) . it is clear that just above the 
first lasing threshold where the spatial hole-burning in 
the inversion is negligible, the threshold lasing mode is 
given by a single CF basis state at the lasing frequency 
satisfying 


??(w) 


7-lA) 

UJ -UJ a + «7_L 


(A9) 


Hence above threshold, the CF basis represents a nu¬ 
merically efficient basis for solving the SALT equations, 
as each lasing mode can be represented by only a few 
CF states. The Sommerfeld condition is implemented by 
matching conditions on the superposition at the last scat¬ 
tering surface to purely outgoing free hermitian solutions. 
The CF equation also facilitates solving the SALT equa¬ 
tions for amplified modes by yielding an additional basis, 
v m (x; uj), j3 m (uj), which are also defined by Eq. (IA 8 I) . but 
satisfy a purely incoming boundary condition instead, 
thus representing a signal with a constant flux into the 
cavity. 

Recently an alternative approach for solving the SALT 
equations has been developed, using a position space ba¬ 
sis, and a truncated simulation region bordered by a per¬ 
fectly matched layer [58|. This follows the conceptual 
approach discussed in detail at the end of Sec. EH The 
Sommerfeld condition is not imposed by matching at a 
last scattering surface, but rather by the presence of the 
PML. In general, the CF methods allow for more elegant 
theoretical analysis, while the position space basis is ulti¬ 
mately expected to be more efficient computationally. In 
both approaches a non-linear iteration scheme will need 
to be implemented to find the correct self-consistent so¬ 
lution in the basis of choice. 

The SALT modal wave equation, Eq. (IA5D is satisfied 
by both lasing modes and injected signals, which inter¬ 
act and compete with each other for gain through the 
non-linear spatial hole-burning denominator of Eq. (IA 6 I) . 
There are two major differences between lasing modes 
and injected signals though; the former must satisfy the 
Sommerfeld radiation condition, while the latter satis¬ 
fies a mixed boundary condition with both incoming and 
outgoing terms. As a direct consequence of this, the fre¬ 
quency of a lasing mode needs to be determined self- 
consistently from the outgoing boundary condition itself, 
while the frequency of an injected signal is externally 
fixed, and the mixed boundary condition can be satisfied 
at any frequency. 


Appendix B: Partially reflecting mirrors in SALT 

Traditionally, SALT has been used to simulate semi¬ 
conductor based microlasers, in which the optical con¬ 
finement for one-dimensional systems is provided by Fres- 
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FIG. 9. Plot of the transmitted intensity of an injected signal 
with frequency Win, into the passive, two-sided cavity with 
partially reflecting mirrors (R = 0.8) on both ends studied 
in Fig. []] Different colors represent different choices of reso¬ 
lution, Ax (red, orange, and green), or the analytic solution 
(blue). 


nel reflection due to the dielectric mismatch across the 
cleaved facets of the device. However, for even the ma¬ 
terials with the largest index of refraction, the optical 
confinement due to Fresnel reflection is weak, R < .36 
for n < 4. Thus to simulate high finesse systems, one 
must either introduce a geometric solution, such as a 
Bragg reflector, or a partially reflecting mirror. In the 
results presented in Sec. 0 we chose to use a partially 
reflecting mirror, as we have the additional motivation of 
simulating a similar system as described by Lang et al. 
and others [3, 153l - l55j . In this appendix we briefly detail 
how to implement a partially reflecting mirror in SALT 
simulations. 

Mathematically, a partially reflecting mirror can be de¬ 
scribed by a delta function in the passive cavity dielectric 
function, 


£c,totai(x) e c (x) T A m 6(x x m ), (Bl) 

m 


where A m represents the strength of the mth mirror at 
x m . For a single mirror, it is straightforward to show 
that the reflection and transmission coefficients are 


R = 



T = 



(B2) 

(B3) 


which demonstrates that the figure of merit for the re¬ 
flectivity of a partially reflecting mirror is fcA. There 
are two different ways to implement this delta-function 
in SALT, either by incorporating the partially reflecting 
mirrors directly into the passive cavity dielectric, as in¬ 
dicated above, or for mirrors at the cavity boundary, ad¬ 
justing the outgoing boundary condition [76| • The former 
method is numerically more stable, and is implemented 
in a discretized system as 



m 


where Sij is the Kronecker delta, x n is the spatial posi¬ 
tion of the nth pixel in the system, and Ax is the res¬ 
olution. Despite the fact that the reflectivity of such a 
mirror is dependent upon £ = kA , when using the CF 
basis to solve the SALT equations, one can choose £ to 
be a constant, and then adjust A accordingly to achieve 
a system with reflectivity independent of the frequency, 
something which is not possible in time-domain simula¬ 
tion techniques. 

We can confirm that the discretized representation of 
a partially reflecting mirror used in conjunction with 
the SALT algorithm is correct by comparing it against 
the analytic result. In Fig. [9] we show the transmission 
resonance for the same passive cavity shown in Fig. [T] 
when calculated analytically (blue), and simulated us¬ 
ing I-SALT with increasing resolution (red, orange, and 
green). The effect of discretizing the system is to shift 
the cavity resonances slightly, and as the resolution is 
increased, the location of the cavity resonance converges 
to the analytic value, as shown. However, regardless of 
the spatial resolution, the Lorentzian shape of the cavity 
resonance remains unchanged. 
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